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1.  Introduction 

This  project  entitled:  'The  structure  of  high-Reynolds-number  turbulent  boundary  layer’ 
was  a  collaborative  effort  between  Yale  (PI:  K.R.  Sreenivasan),  Princeton  (PI:  A.J.  Smits)  and 
Penn.  State  (PI:  J.  Brasseur).  The  main  point  of  the  proposal  was  that  the  structure  of  the  high- 
Reynolds-number  turbulent  boundary  layer  may  differ  substantially  from  that  at  moderate 
Reynolds  numbers,  and  is  therefore  worth  a  special  look. 

The  emphasis  in  this  project  -  which  involved  experiment,  simulations  and  data  analysis  - 
was  to  be  on  aspects  such  as:  (a)  in  what  respects  does  the  boundary  layer  at  high  Reynolds 
numbers  differ  from  that  at  low  Reynolds  numbers?  (b)  What  mechanisms  in  particular  contribute 
to  turbulent  transport  of  heat,  mass  and  momentum  in  the  logarithmic  region?  (c)  What  sets  the 
scales  for  transport  in  the  boundary  layer?  (d)  Is  the  transport  primarily  in  the  form  of  some  small 
number  of  very  active  events  occupying  small  amounts  of  time  (and  space),  or  does  it  occur  over 
most  of  the  time  (and  occupy  most  of  the  space)  in  a  fashion  that  is  not  spectacular?  (e)  How  can 
the  transport  phenomenon  be  modelled  for  use  in  sub-grid  scale  models?  (0  What  analogies  from 
other  branches  of  physical  sciences  can  be  brought  to  bear  on  the  transport  issue?  Although 
questions  concerning  energy  dissipation,  velocity  increments,  zero-crossings,  fractal  dimensions, 
were  to  be  addressed,  the  central  part  of  the  work  was  to  be  on  turbulent  transport.  Progress  made 
at  Yale  on  these  and  related  questions  will  be  described  here  briefly. 


2.  Notes  on  the  Yale  work 

(a)  Data  were  acquired  in  the  laboratory  as  well  as  atmospheric  boundary  layer. 
Atmospheric  data  were  obtained  at  a  few  selected  heights  above  a  wheat-field  canopy.  They  were 
primarily  in  the  form  of  'single  point'  measurements.  The  quantities  measured  were  velocity 
components  (two),  the  Reynolds  shear  stress,  vorticity  (streamwise  component)  and  temperature. 
'Flow  visualization'  was  done  by  using  smoke  bombs  (generally  lasting  for  about  three  minutes)  at 
several  heights.  Laboratory  data  were  comparable. 
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(b)  The  multiplier  method  was  develeoped  for  displaying  scale-similarity  in  the  momentum 
transport  process  at  high  Reynolds  numbers;  new  schemes  for  the  fractal  analysis  of  turbulent 
velocity  and  temperature  data  in  the  atmospheric  surface  layer  were  also  developed  and  used. 

(c)  The  issue  of  incomplete  similarity  (Barenblatt,  J.  Fluid  Mech.  248,  513,  1993)  was 
examined  in  so  far  it  bears  on  the  mean  velocity  distribution  in  high-Reynolds-number  pipe  flows 
and  boundary  layers. 

(d)  The  Reynolds  number  effects  on  the  inner-outer  interactions  were  investigated, 
especially  on  the  Reynolds  shear  stress. 

(e)  Several  statistical  properties  were  examined  over  a  wide  range  of  Reynolds  numbers; 
these  included  the  probability  density  function  (pdf)  of  the  zero-crossing  intervals  in  turbulent 
signals,  and  the  pdf  of  velocity  increments. 

(f)  The  scalar  field  as  well  as  velocity  data  were  mapped  out  in  two-dimensional  planes  (the 
latter  using  PIV  techniques)  in  a  moderate-Reynolds-number  boundary  layer.  The  reason  for  the 
moderate-Reynolds-number  measurements  is  that  it  can  be  probed  in  much  greater  detail  than  is 
possible  (for  technical  reasons)  at  high  Reynolds  numbers,  and  that  the  data  so  obtained  may  serve 
as  a  basis  for  comparison  with  the  high-Reynolds-number  data.  It  would  similarly  serve  well  for 
the  examination  of  the  Reynolds  number  effects. 

3.  Some  further  details 

(a)  Investigation  of  the  skin  friction  data  in  pipes  and  boundary  layers  indicated  that  there 
might  be  a  reasonably  well-defined  'transition'  between  the  low  Reynolds  number  behavior  and  the 
high  Reynolds  number  behavior.  Velocity  profiles  as  well  as  pressure  tluctuations  seemed  to 
support  the  notion  that  a  weak,  or  second-order,  transition  may  indeed  occur  at  a  fairly  high 
Reynolds  number.  This  Reynolds  number  for  the  boundary  layer  is  of  the  order  of  6000  in  terms 
of  the  momentum  thickness,  or  for  the  pipe  of  the  order  of  about  105  in  terms  of  the  pipe  diameter. 
The  following  physical  effects  were  considered  as  possible  causes  for  this  'transition',  and  their 
relative  effects  were  investigated. 


(i)  It  is  possible  that  these  effects  could  arise  because  the  importance  of  the  viscous- 
influenced  interface  in  the  outer  region  could  become  weak  beyond  some  Reynolds  number. 
Fractal  estimates  and  zero-crossing  measurements  precluded  such  a  possibility. 

(ii)  It  is  possible  that  the  localization  of  the  energy  dissipation  in  the  outer  region  may  be  the 
cause.  But  multiffactal  estimates  show  that  this  is  an  unlikely  because  the  localization  is  extremely 
slow  in  Reynolds  number. 

(iii)  It  is  possible  that  these  effects  could  come  from  the  altered  circumstances  of  the 
inner/outer  interaction.  It  was  first  shown  by  an  analysis  of  rough-pipe  data  that  the  really 
important  parameter  is  not  the  Reynolds  number  but  the  ratio  of  the  inner  scale  to  the  outer  scale. 
This  ratio  is  the  same  as  the  Reynolds  number  in  the  smooth-wall  flows,  but  it  can  be  prescribed 
independent  of  the  Reynolds  number  for  flows  with  rough  walls.  Thus,  the  so-called  Reynolds 
number  effects  are  really  the  effects  of  the  outer/inner  scale  ratio.  This  also  suggests  that  a  study  of 
the  high-Reynolds-number  dynamics  is  not  ideally  studied  in  the  atmospheric  boundary  layer  over 
rough  terrain  where,  even  though  the  Reynolds  number  is  la-ge,  the  effective  scale  ratio  is  not  so 
large. 

(b)  It  is  shown  that  for  large  values  of  the  outer/inner  scale  ratio,  as  opposed  to  the  case  of 
moderate  values  of  this  ratio,  the  velocity-conditioned  Reynolds  shear  stress  changes  its  character. 

(c)  It  is  shown  that  multiplicative  representations  of  the  Reynolds  stress,  in  spite  of  some 
advantages  of  interpretation,  has  several  problems  in  the  logarithmic  region  of  the  boundary  layer. 
In  any  case,  scale-similar  ranges  for  the  multiplier  distribution  of  the  Reynolds  shear  stress  have 
been  identified. 

(d)  It  is  shown  by  investigating  the  instantaneous  momentum  fluxes  that,  on  the  average, 
the  contributions  to  the  Reynolds  stress  comes  from  medium  amplitude  velocity  fluctuations  of  the 
order  of  1.5  root-mean -square  value.  Thus  very  large  velocity  amplitudes  do  not  contribute  much 
to  the  Reynolds  stress. 


4.  Publications 

Most  of  this  work  can  be  found  in  the  following  references: 

(a)  P.  Kailasnath,  K.R.  Sreenivasan  and  G.  Stolovitzky,  "Probability  density  of  velocity 
increments  in  turbulent  flows",  Phys.  Rev.  Lett.  68,  2766,  1992 

(b)  A.B.  Chhabra  and  K.R.  Sreenivasan,  "Scale-invariant  multiplier  distributions  in 
turbulence",  Phys.  Rev.  Lett.  68,  2762,  1992 

(c)  K.R.  Sreenivasan  and  A.  Juneja,  "Fractal  dimensions  of  time  series  in  turbulent  flows", 
submitted  to  J.  Fluid  Mech.,  1992 

(d)  P.  Kailasnath  and  K.R.  Sreenivasan,  "Zero-crossings  of  velocity  fluctuations  in 
turbulent  boundary  layers",  submitted  to  Phys.  Fluids  A.,  1992 

(e)  P.  Kailasnath,  "Studies  of  Reynolds  number  effects  in  turbulent  boundary  layers", 
Ph.D.  thesis,  Yale  University,  June  1993 

The  papers  (a)  and  (b)  are  enclosed.  Items  (c),  (d)  and  (e)  are  in  the  process  of  being 
published  and,  once  published,  copies  will  be  sent  to  AFOSR.  Most  of  the  material  central  to  the 
proposal  is  contained  in  item  (e). 

5.  Collaboration  with  Princeton  and  Penn.  State 

i 

The  collaborative  aspect  of  the  work  proposed  between  Yale,  Penn.  State  and  Princeton. 
This  collaboration  is  still  continuing  (even  though  the  project  has  formally  ended),  and  it  is 
expected  that  a  few  significant  publications  will  arise  from  it.  The  collaboration  has  taken  the  form 
of  exchange  of  data,  mutual  visits  and  coordination  of  flow  conditions  between  simultaions  and 
experiments.  In  terms  of  concrete  results,  we  may  cite  the  following  aspects: 

(a)  Comparison  of  the  inertial  region  of  a  high-Reynolds-number  turbulent  boundary  layer 
with  homogeneous  shear  flow; 

(b)  Wavelet  analysis  of  the  experimental  as  well  as  numerically  generated  data  in  boundary 

layers; 


(c)  Comparison  between  the  simulations  and  experiments  in  moderate-Reynolds-number 
boundary  layers. 

A  draft  of  a  paper  on  item  (a)  is  available,  and  (a)  as  well  as  (b)  have  been  presented  at 
APS/DFD  meetings,  as  well  as  elsewhere.  It  will  be  several  more  months  before  journal 
publications  materialize. 

6.  Aspects  of  the  work  that  were  not  successful 

(a)  We  had  proposed  to  make  measurements  in  high-Reynolds-number  boundary  layer  in 
the  new  wind-tunnel  facility  that  was  to  have  been  completed  at  Illinois  Institute  of  Technology 
about  two  or  three  years  ago.  Unfortunately,  this  was  not  possible  because  the  wind  tunnel  is  still 
unfinished  as  of  this  writing.  So  we  made  progress  on  all  other  fronts  available  to  us:  acquisition 
and  analysis  of  data  in  the  high-Reynolds-number  atmospheric  surface  layer  (already  mentioned), 
examination  of  the  high-Reynolds-number  mean  velocity  data  acquired  at  APL  of  Penn.  State,  use 
of  some  relatively  high-Reynolds-number  data  acquired  at  Princeton,  and  a  study  of  the  moderate 
Reynolds  number  data  acquired  by  us.  To  supplement  the  analysis,  we  also  examined  all  previous 
data  obtained  elsewhere,  in  both  boundary  layers  and  pipe  flows. 

(b)  We  had  planned  on  obtaining  the  scalar  and  velocity  fields  simultaneously  at  moderate 
Reynolds  numbers,  but  did  not  have  the  time  to  complete  these  measurements  before  the  project 
came  to  an  end,  and  the  post-doctoral  fellow  working  on  the  problem  had  to  be  terminated. 
Unfortunately,  much  time  and  effort  was  invested  in  this  work.  The  consolation  is  that  the 
expertise  acquired  will  be  of  use  in  subsequent  work. 

7.  Personnel 

Besides  the  PI,  much  of  the  work  on  this  project  was  done  by  a  post-doctoral  fellow.  Dr. 
A.K.  Suri  (who  made  most  of  the  scalar  and  PIV  measurements  in  moderate-Reynolds-number 
boundary  layers),  and  a  graduate  student,  Mr.  P.  Kailasnath  (who  took  a  large  share  of  the 
responsibility  for  single-point  measurements  in  laboratory  as  well  as  atmospheric  boundary 
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layers).  Kailasnath  has  just  completed  his  Ph.D.  thesis  and  has  taken  up  a  position  at  the  Yale 
medical  School.  Off  and  on,  as  demanded  by  circumstances,  other  people  were  employed  on 
short-term  basis.  Some  of  the  work  related  to  the  vorticity  measurements  in  the  atmospheric 
boundary  layer  are  contained  in  the  Ph.D.  dissertation  of  Mark  S.  Fan.  Mark's  thesis  was  only 
partly  supported  by  another  grant  AFOSR-87-Ol  16.  He  now  works  for  NASA. 

A  final  note:  There  are  several  aspects  of  this  work  that  are  yet  to  be  digested  fully,  and  it 
is  our  responsible  belief  that  more  new  results  will  appear  out  of  this  work  in  due  course. 
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Scale-Invariant  Multiplier  Distributions  in  Turbulence 

Ashvin  B.  Chhabra<l)(2)  and  K.  R.  Sreenivasan(l> 

01  Mason  Laboratory.  Yale  University.  New  Haven.  Connecticut  06520-2159 
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A  family  of  scale-invariant,  base-dependent,  multiplier  distributions  is  measured  for  the  turbulence 
dissipation  field  in  the  atmospheric  surface  layer.  The  existence  of  these  distributions  implies  the  ex¬ 
istence  of  the  more  traditional  multifractal  scaling  functions,  and  we  compute  both  positive  and  negative 
parts  of  the  /(a)  curve.  The  results  support  the  conjecture  of  universality  in  the  scaling  properties  of 
small-scale  turbulence.  A  simple  cascade  model  based  on  the  measured  multiplier  distributions  is  shown 
to  possess  several  advantages  over  previously  considered  models. 

PACS  numbers:  47.25.-c.  02.50. +s.  03.40.Gc.  05.45.+b 


The  multifractal  formalism  [ I  —5]  le.g.,  the  /(a)  func¬ 
tion!  compactly  describes  the  scaling  properties  of  mea¬ 
sures  which  arise  in  a  variety  of  problems  such  as  chaotic 
dynamical  systems  [6],  diffusion-limited  aggregation 
(DLA)  17],  and  dissipationlike  quantities  in  fully  devel¬ 
oped  turbulence  [8].  The  /(a)  function  is  a  macroscopic 
feature  of  such  measures,  and  a  quest  for  deeper  under¬ 
standing  of  the  underlying  physics  has  led  to  microscopic 
descriptions.  A  successful  example  is  the  Feigenbaum 
scaling  function  (FSF)  19]  appropriate  to  the  onset  of 
chaos  in  the  period-doubling  route.  Basic  to  the  FSF  are 
the  contraction  ratios  (or  multipliers),  which  describe 
how  distances  between  the  nearest-neighbor  iterates  scale 
with  increasing  levels  of  refinement.  The  FSF  organizes 
these  multipliers  according  to  the  natural  time  of  the  sys¬ 
tem  described  by  the  closest  return  times,  and  compactly 
describes  the  local  scaling.  The  FSF  description  is  much 
richer  than  that  embodied  in  /(a),  and  the  latter  is  easily 
computable  from  the  former.  Further,  the  /(a)  descrip¬ 
tion  is  degenerate  and  a  variety  of  scaling  functions  lead 
to  the  same  /(a)  function  1 10,1 1]. 

Many  muitifractai  measures,  sucl.  as  the  spatial  distri¬ 
bution  of  the  energy  dissipation  rate  in  fully  developed 
turbulence,  or  the  harmonic  measure  of  the  DLA  struc¬ 
tures,  display  statistical  properties  that  are  different  from 
those  of  deterministic  systems.  Two  such  properties  are 
the  sample-to-sample  fluctuations  of  the  /(a)  function, 
and  the  existence  of  negative  dimensions —  both  of  which 
represent  the  same  underlying  phenomenon  [12,13].  Re¬ 
call  that,  if  P,(e)  is  the  measure  in  the  ith  box  of  size  e, 
one  can  decompose  the  multifractal  measure  into  in¬ 
terwoven  sets  of  varying  singularity  strengths  a  (where 
I/*, (c) whose  fractal  dimensions  /(a)  are  de¬ 
fined  by  N(a)  being  the  number  of  singu¬ 

larities  of  strength  a.  Figure  I  shows  how  /(q**  2)  and 
a(q~2),  computed  from  the  turbulence  dissipation  data 
using  the  direct  method  of  Ref.  (14],  vary  from  one  sam¬ 
ple  to  another;  each  sample  is  approximately  one  macro¬ 
scale  in  extent  and  produces  extended  and  unambiguous 
scaling.  The  sample-to-sample  fluctuations  (I Si  are 
much  larger  than  the  least-squares  errors  in  calculating 


f(q)  and  a(q)  from  individual  samples.  In  addition,  for 
any  given  q,  f(q)  and  a(q)  from  different  samples  are 
correlated  and  fall  along  a  thin  band  [16]. 

These  phenomena  are  absent  in  deterministic  processes, 
and  their  occurrence  reflects  an  inherently  probabilistic 
dynamics.  For  such  cases,  a  complete  specification  of  the 
scaling  properties,  even  at  the  level  of  /(a),  requires 
measuring  both  positive  and  negative  parts  of  f(a) 
[12,13,17-20].  Negative  f(a)  for  turbulent  flows  at 
moderate  Reynolds  numbers  have  already  been  obtained 
in  Ref.  [20]  but,  for  atmospheric  flows,  their  computation 
by  conventional  box-counting  methods  requires  an  enor¬ 
mous  amount  of  data,  involving  perhaps  several  years  of 
data  acquisition  [12,13,20].  Recently  [12,13,21],  a 


FIG.  I.  Sample-to-sample  fluctuations  of  f(q“ 2)  vs 
a(q  — 2)  for  ten  samples  (7200  points  each,  roughly  equivalent 
to  an  integral  scale)  from  an  atmospheric  boundary  layer 
(squares).  The  circles  are  from  a  simulation  of  a  binary  cas¬ 
cade  model  with  multipliers  chosen  randomly  from  a  triangular 
distribution.  This  model  is  described  later  in  the  text  (see  also 
Fig.  2). 


2762 


€>  1992  The  American  Physical  Society 


Volume  68,  Number  18 


PHYSICAL  REVIEW  LETTERS 


4  May  1991a 


method  based  on  scale-invariant  multiplier  distributions 
has  been  proposed  as  being  exponentially  faster  and  more 
accurate  than  box-counting  methods.  Utilizing  this 
method,  we  were  able  to  compute  negative  dimensions  for 
atmospheric  flows  as  well  as  a  variety  of  model  cascade 
processes. 

This  paper  has  several  related  objectives.  As  a  con¬ 
tinuation  of  the  notion  that  scale-invariant  multiplier  dis¬ 
tributions  are  fundamental  functions,  we  compute  them 
for  atmospheric  turbulence  by  assuming  that  the  observed 
scaling  could  result  from  cascades  of  any  base  (binary, 
ternary,  etc.).  Using  these  distributions,  we  construct 
several  base-independent  functions  such  as  f(a).  We  then 
show  that  the  f(a)  functions  computed  from  different 
multiplier  distributions  agree  well  with  each  other,  as  well 
as  with  those  computed  in  moderate-Reynolds-number 
flows  by  box-counting  methods.  Finally,  we  propose  a 
relatively  simple  cascade  model  that  reproduces  not  only 
the  observed  positive  and  negative  parts  of  /(a)  and 
sample-to-sample  fluctuations,  but  also  the  stretched  ex¬ 
ponential  tails  in  the  probability  distribution  of  tur¬ 
bulence  dissipation.  The  superiority  of  the  multiplier 
method  allows  us  to  compute  /(a)  with  greater  accuracy 
and  to  assess,  for  the  first  time,  the  universality  conjec¬ 
ture  about  scaling  properties  in  the  negative-dimension 
range. 

To  compute  the  scale-invariant  multiplier  distributions 
one  first  construes  a  measure  which,  in  the  present  case, 
is  the  energy  dissipation  rate  represented  by  the  square  of 
the  derivative  of  a  component  of  the  velocity.  The  mea¬ 
sure  is  then  covered  by  boxes  of  uniform  size.  Each  of 
these  boxes  is  then  subdivided  into  a  number  of  boxes,  a , 
and  the  ratios  of  the  measures  in  the  original  box  to  those 
in  the  smaller  sub-boxes  are  computed.  A  histogram  of 
these  ratios  is  then  Pa(Ma)  Figure  2  shows  Pa(Ma)  ob¬ 
tained  by  assuming  cascades  of  bases  a  —2,  3,  and  5.  The 
larger  symbols  show  an  average  over  steps  involving  com¬ 
parisons  between  boxes  of  size  m  and  those  of  size  ma, 
where  m  ranged  from  50  to  1000.  The  shape  of  the  dis¬ 
tribution  remains  invariant  for  the  inertial  range  of 
scales.  (For  the  smallest  scales,  the  distributions  have  a 
concave  shape.  This  concavity  is  related  to  the  diver¬ 
gence  of  moments  [22]  and  will  be  discussed  elsewhere. 
For  very  large  boxes,  multiplier  distributions  become 
flatter,  as  would  be  the  case  for  random  measures.) 

The  multiplier  distributions  P„(Ma )  are  more  basic 
than  the  /(«)  function,  and  were  introduced  in  the  con¬ 
text  of  turbulence  by  Novikov  111,  and  measured  by  Van 
Atta  and  Yeh  [23].  However,  their  importance  with  re¬ 
gard  to  multifractals  in  turbulence  has  subsequently  been 
ignored.  One  can  compute  from  the  multiplier  distribu¬ 
tions  both  the  positive  and  negative  parts  of  the  /(a) 
function,  i.e.,  they  contain  information  on  the  asymptotic 
scaling  properties  of  a  measure  and  on  fluctuations  of 
scaling  properties  for  samples  of  finite  size.  In  addition, 
even  in  instances  where  high-order  moments  diverge 
[2,22],  P„(Ma)  remains  well  defined.  Finally,  the  /(«) 


FIG.  2.  Multiplier  distributions  P.(Ma )  for  bases  (from 
right  to  left)  a  “2,  3,  and  5.  The  larger  symbols  show  averages 
over  steps  comparing  measures  in  boxes  of  size  m  with  those  in 
size  ma,  SO  <  m  <  1000  in  units  of  the  Kolmogorov  scale.  The 
smaller  symbOi„  show  the  distributions  obtained  for  different 
values  of  m  (50,  80.  150.  200.  400.  and  1000)  The  solid  line  is 
a  useful  approximation  to  the  a  “2  distribution 


function  may  extend  over  (  —  oo.oo)  whereas  Pa(Ma )  is  a 
compact  function  defined  on  Ma  €  [0,1].  The  base 
dependency  of  P„(Ma)  can  be  scaled  out  because  the 
multiplier  distributions  corresponding  to  different  bases 
are  related  by  convolution,  provided  the  multipliers  at 
successive  cascades  are  uncorrelated.  With  this  assump¬ 
tion,  several  base-independent  functions  can  be  construct¬ 
ed  from  the  multiplier  distributions. 

Consider  as  an  example  the  Mellin  transforms  M  of 
Pa(.M0)  and  of  its  convolution  PiOV),  where  A  is  the 
product  of  two  multipliers  picked  according  to  their  prob¬ 
abilities  Pa(Ma).  Since  the  Mellin  transform  of  a  convo¬ 
lution  is  the  product  of  the  individual  Mellin  transforms, 
we  have 

M \Pt-a (Mh )}  - M  \Pa (M„ ))  M \Pa(Ma )! 

-\M\Pa(M.)W  (I) 

It  is  clear  that,  in  general,  the  exponent  in  the  last  term 
of  Eq.  (I)  is  simply  log(b)/log(a),  the  number  of  times 
the  variable  is  being  convolved.  So  far,  for  any  two 
different  bases  we  have 

[M \P„(Ma )i ]  ,,/lo*<a,!  - [M {Pb(Mh )}  1 ::  ,~*(c  'l  (2) 


On  evaluating  the  Mellin  transform  and  taking  loga¬ 
rithms  of  both  sides,  we  get  [I] 


log ((*/„)«)  _  log {(Mt)') 
log(a)  log(b) 


[r(^)  +  f)0] 


(3) 


The  other  scaling  functions f(q),  a(q),  and  Dq  are  simply 
related  to  r(q),  and  so  they  too  can  be  easily  derived 
from  the  multiplier  distributions.  Further,  one  can  derive 
equivalent  scaling  functions  by  using  Laplace  or  Fourier 
transforms. 

We  now  compute  the  f(a)  function  from  the  multiplier 
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FIG.  3.  The /(a)  function  from  the  atmospheric  data  for  FIG.  4.  The  f(a)  function  for  the  atmospheric  data  for  bases 

bases  2.  3.  4,  5,  6.  8.  and  10,  after  the  base-dependent  prefactor  2  and  10.  Also  reproduced  (from  Ref.  (20])  are  the  data  from 

has  been  removed.  The  good  collapse  implies  that  the  assump-  a  variety  of  laboratory  flows  (larger  symbols).  The  small  dia- 
tion  of  random  multiplicative  processes  is  essentially  correct.  monds  are  from  compulations  l28l  using  multipliers  from  the 

and  that  no  single  base  is  preferred.  triangular  distribution  shown  in  Fig.  2.  The  agreement  between 

/(a)  functions  for  atmospheric  turbulence  and  for  laboratory 
flows  supports  the  notion  of  universality  of  scaling  properties  at 
distributions  of  Fig.  2,  and  verify  that  the  multifractal  different  Reynolds  numbers  even  of  the  most  rare  events, 
scaling  functions  derived  from  different-base  multiplier 
distributions  indeed  collapse.  While  doing  so.  it  is  impor¬ 
tant  to  account  (  for  ^the  existence  of  a  prefactor  in  parisoninthenegative/(a)rcgionshowsexcellentagree- 
{(Ma)q) ~C(q)a  '  ’  From  Eq.  (3),  we  have  ment  between  the  laboratory  and  atmospheric  scaling 

t „(q)  +  Do~  t  (q) -log[C(q)l/log(a)  +  Do.  but  the  pre-  properties,  supporting  the  conjecture  of  universal  scaling 

factor  can  be  eliminated  by  computing  ra(q)  for  two  for  even  the  rare  events  in  the  small-scale  velocity  field  in 

different  bases.  One  can  then  obtain  the  correct  base-  fully  developed  turbulence. 

dependent  exponents  for  an  arbitrary  base.  Internal  con-  In  the  past  few  years  a  variety  of  cascade  models  such 
sistency  requires  that  all  of  them  should  collapse  onto  a  as  the  log-normal  model  [24].  beta  model  l25l,  random 

single  curve.  Figure  3  shows  just  such  a  collapse  for  beta  model  [26),  and  p  model  18)  have  been  proposed  to 

bases  a  —2,  3,  4,  5,  6,  8,  and  10  (where  the  prefactor  was  mimic  intermittency  in  turbulence.  None  of  these  models 

computed  using  bases  2  and  4).  The  collapse  strongly  in-  displays  sample-to-sample  fluctuations  or  negative  dimen- 

dicates  that  there  is  no  preferred  base  as  far  as  it  con-  sions;  nor  do  they  explain  the  multiplier  distributions 

cents  the  scaling  properties.  shown  in  Fig.  2.  One  can.  however,  construct  an  entire 

Finally,  Fig.  4  compares  the  scaling  exponents  for  at-  family  of  simple  cascade  models  which  display  all  of  the 

mospheric  flows  with  those  computed  [20l  from  laborato-  above  features  and  whose  /(a)  functions  agree  well  with 

ry  data.  The  latter  have  been  computed  using  conven-  experiment.  The  simplest  model  would  be  a  binary  model 

tional  box-counting  methods  and  thus  do  not  assume  un-  with  multipliers  picked  randomlv  from  a  triangular  distri- 

correlated  multiplicative  processes.  The  exponents  for  at-  bution  (shown  by  the  solid  line  in  Fig.  2)  which  is  a  rough 

mospheric  flows  have  been  computed  using  the  multiplier  approximation  to  the  mean  distribution  for  a  binary  cas- 

method  which  uses  just  such  an  assumption.  The  excel-  cade.  Figure  4  shows  that  the  /(a)  function  for  this 

lent  agreement  with  conventional  box  counting  (in  the  model  is  in  good  agreement  with  other  experimental  data, 

range  where  the  latter  is  capable  of  yielding  exponents)  As  already  remarked,  this  cascade  model  displays  the 

indicates  the  correctness  of  our  convolution  arguments  right  behavior  with  respect  to  sample-to-sample  fluctua- 

bascd  on  uncorrelated  multiplicative  processes.  (Such  an  tions  shown  by  circles  in  Fig.  I .  The  model  also  repro¬ 
approach  would  fail  for  the  period-doubling  attractor  be-  duces  the  observed  [20,27)  stretched  exponential  tails, 

cause  the  contraction  ratios  in  that  case  are  highly  corre-  P(€)~cxpf-0(e)l/Ji,  in  the  probability  distributions  of 

lated,  thus  invalidating  the  simple  convolution  arguments  the  energy  dissipation  rate  e.  However,  it  cannot  address 

used  here.)  We  thus  observe  a  nontrivial  and  interesting  issues  such  as  the  divergence  of  high-order  moments, 

empirical  fact  that  the  successive  multipliers  in  cascades.  We  thank  J.  Eggers,  C.  J.  Evertsz.  R.  V.  Jensen,  P.  W. 
which  give  rise  to  the  observed  intermittency  in  tur-  Jones,  L.  P.  Kadanoff.  B.  B  Mandelbrot,  C  Mcneveau. 

bulencc.  seem  to  be  essentially  uncorrelated.  The  com-  M.  Nelkin,  N.  Read,  and  G.  Stolovitzky  for  useful  dis- 
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Measurements  have  been  made  of  the  probability  density  function  (PDF)  of  velocity  increments  A u(r) 
for  a  wide  range  of  separation  distances  r.  Stretched  exponentials  provide  good  working  approximations 
to  the  tails  of  the  PDF.  The  stretching  exponent  varies  monotonically  from  O.S  for  r  in  the  dissipation 
range  to  2  for  r  in  the  integral  scale  range.  Theoretical  forms  based  on  multifractal  notions  of  tur¬ 
bulence  agree  well  with  the  measured  PDFs.  When  the  largest  scales  in  the  velocity  u  are  filtered  out, 
the  PDF  of  A u(r)  becomes  symmetric  and,  for  large  r,  close  to  exponential. 


PACS  numbers  47.25. -c.  02.50. +s.  03.40.Gc.  05.45.+b 

Much  work  in  the  turbulence  literature  (1-4]  has  been 
devoted  to  the  determination  of  the  scaling  properties  of 
the  structure  functions  ([Au(/’)1"),  where  n  is  a  positive 
integer  and  Au(r)  is  the  velocity  increment  between  two 
spatial  locations  which  are  a  distance  r  apart.  It  has  re¬ 
cently  (51  been  emphasized  that  a  better  strategy  may  be 
to  focus  on  the  probability  density  functions  (PDF  )  of 
Ah(z),  p^(Nu),  rather  than  on  the  collection  of  mo¬ 
ments.  Accordingly,  this  Letter  is  concerned  with  the 
PDFs  of  Au(r)  and  of  velocit"  derivatives,  and  has  three 
purposes.  First,  it  provides  experimental  data  on 
/^(Aufr))  for  a  large  range  of  separation  distances  r 
spanning  the  dissipation  range  on  one  end  and  the  in¬ 
tegral  scale  range  on  the  other.  The  data  show  that  the 
PDFs  are  square-root  exponential  for  r  in  the  dissipation 
range  and  Gaussian  for  r  on  the  order  of  the  correlation 
(or  the  integral)  scale.  Second,  it  provides  for  the  PDFs 
a  theoretical  expression  based  on  the  multifractal  picture 
for  turbulence  dynamics.  The  third  aspect  is  related  to 
the  exact  result  known  for  velocity  increments  [6],  name¬ 
ly,  that  in  locally  homogeneous  turbulence,  the  third- 
order  structure  function  in  the  inertial  range  obeys  the  re¬ 
lation 

<[A«(r)]3>--f<e>r.  (I) 

While  the  implications  of  the  nonzero  value  of  ((Au(r)]3) 
have  not  been  understood  fully,  it  has  been  shown  recent¬ 
ly  (7]  that  the  removal  of  the  largest  scales  by  high-pass 
filtering  renders  the  PDF  of  Ah  symmetric.  This  Letter 
provides  experimental  results  on  the  p &  when  the  low- 
frequency  (or  large  scale)  components  of  u  have  been 
filtered  out.  For  complementary  results  on  PDFs  of  ve¬ 
locity  derivatives  and  increments,  see  Refs.  (8]  and  (9], 
We  shall  subsequently  return  to  some  of  this  work. 

Measurements  were  made  in  the  atmospheric  surface 
layer  about  6  m  above  a  wheat  canopy  in  the  Connecticut 
Agricultural  Research  Station.  Data  were  also  acquired 
about  2  m  above  the  roof  of  a  four-story  building.  The 
laboratory  data  were  acquired  at  a  height  of  0.2 S,  where 
5  is  the  thickness  of  the  boundary  layer,  over  a  smooth 
flat  plate.  The  boundary-layer-thickness  Reynolds  num¬ 
ber  at  the  measuring  station  was  32000.  Velocity  fluc¬ 
tuations  were  measured  using  the  standard  hot-wire  (S 
pm  diam,  0.6  mm  length)  velocimeter  operated  in  the 


constant-temperature  mode  on  a  DISA  S5M0I  anemom¬ 
eter.  The  anemometer  voltage  was  digitized  on  a  1 2-bit 
digitizer  and  linearized  before  further  processing.  Veloci¬ 
ty  derivatives  were  obtained  by  central  differencing  of  the 
data.  Taylor’s  frozen  Row  hypothesis  was  used  in  inter¬ 
preting  time  intervals  as  space  intervals.  The  precise  lim¬ 
itations  of  this  hypothesis  are  unclear  (in  spite  of  much 
work),  especially  for  the  tails  of  the  PDF,  but  it  should  be 
noted  that  the  mean  convection  velocity  in  the  present  ex¬ 
periments  was  about  15  times  larger  than  the  standard 
deviation  of  the  fluctuating  velocity.  For  obtaining  the 
PDF  of  the  filtered  data,  a  linear  phase  filter  with  excel¬ 
lent  cutoff  characteristics  (10]  was  used. 

We  have  fitted  stretched  exponentials,  ptu, 
~-exp(  —  alAnl"),  to  the  tails  of  p^fAu).  The  inset  to 
Fig.  I  demonstrates,  for  one  arbitrarily  chosen  r,  that 


FIG.  I.  The  experimentally  determined  stretching  exponent 
m  in  p*.(Aw(r))~ expl —  a(Au)"],  plotted  as  a  function  of 
r/A.  O  and  □  are  for  two  different  sets  of  atmospheric  data.  The 
Taylor  microscale  Reynolds  number  R»  was  on  the  order  of 
1500  for  both.  The  laboratory  data,  not  plotted  here,  show  a 
similar  trend  and  roughly  coincide  with  the  atmospheric  data. 
Au(r)  is  obtained  by  taking  the  velocity  differences  separated 
by  a  time  difference  A t,  and  interpreting  r  “  —  AtU.  where  U  is 
(he  mean  velocity  of  the  flow  at  the  measurement  station 
(Taylor's  frozen  flow  hypothesis).  The  integral  scale  A  was 
determined  by  obtaining  the  area  under  the  autocorrelation 
function  of  u.  and  converting  it  to  a  length  scale  by  Taylor’s  hy¬ 
pothesis.  Inset:  The  stretched  exponential  is  a  good  approxima¬ 
tion  to  typical  experimental  data  for  both  sides  of  the  distribu¬ 
tion.  In  this  figure  and  others  in  this  paper,  A u  is  normalized  by 
its  root-mean-square  value. 
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there  is  an  extensive  region  of  the  PDF  to  which  a 
stretched-exponential  fit  is  good.  It  follows  from  Eq.  (I) 
that  there  must  be  a  certain  asymmetry  between  the  two 
tails  of  the  distribution,  but  this  asymmetry  is  not  very 
large.  In  fact.  Fig.  I  shows  that  the  differences  between 
the  two  tails  of  the  distribution,  insofar  as  they  relate  to 
the  stretched-exponential  fits,  are  small.  We  shall  there¬ 
fore  momentarily  ignore  this  asymmetry  and  return  to  it 
later.  The  empirically  determined  stretching  exponents 
are  plotted  in  Fig.  1  as  a  function  of  r  for  two  sets  of  at¬ 
mospheric  data.  While  square-root-exponential  fits  are 
good  for  r  in  the  dissipation  range,  Gaussian  fits  are  ap¬ 
propriate  for  r— Z-o-lOA,  where  A  is  the  autocorrelation 
length  scale  of  u  (see  caption  to  Fig.  I).  This  latter  result 
is  not  surprising  because  velocities  at  two  widely  separat¬ 
ed  points  become  independent  of  each  other.  The  ex¬ 
ponent  m  increases  monotonically  through  the  inertial 
range.  Laboratory  data,  not  displayed  here,  show  a  simi¬ 
lar  behavior. 

It  should  be  noted  that  a  somewhat  similar  effort  for 


temperature  increments  in  a  Rayleigh-Benard  convection 
experiment  has  been  made  in  Ref.  [I  I]. 

In  a  recent  paper  [9],  Benzi  et  al.  derived  an  expression 
for  the  PDF  of  velocity  increments.  They  assumed  that 
A u(r),  for  r  in  the  inertial  range,  is  given  by  the  random 
P  model  [12], 

Au(r)*“AuoWL)l/3n^rl/3 .  (2) 

i-i 

where  Auo  is  the  characteristic  velocity  increment  on  the 
macroscopic  length  Lo ,  and  the  Pi's  are  identically  distri¬ 
buted  independent  random  variables.  Benzi  et  al.  [9] 
used  a  special  case  in  which  the  probability  density  of  P 
was  assumed  to  be  given  by 

p,(p)-aS(\-p)  +  (l-a)S(B~p) ,  (3) 

with  a“  i  and  B“j,  in  conformity  with  experiments 
[2,3],  Under  the  further  assumption  that  Amo  is  normally 
.  distributed  (see  Fig.  I ),  Benzi  et  al.  [9]  showed  that  the 
’  PDF  of  Am  is  given  by 


^(Aufr))-  Z  Cfa"-‘(1  -a)kBk— -1-^ 
*-o  (2  xal^r'1 


exp 


<*k. 


<A  ufruHr/L)m 
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+  ll-[a  +  (|-a)«]"U(A«), 
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Ub) 


Here,  n  is  the  number  of  steps  assumed  to  occur  in  a  cas¬ 
cade  before  reaching  a  scale  r  and  is  given,  for  a  binary 
cascade,  by  n~\of,i(Lo/r).  Notice  that  the  second  term 
on  the  right-hand  side  of  Eq.  (4a)  accounts  for  inactive 
eddies,  i.e.,  spatial  regions  where  the  dissipation  of  energy 
is  zero.  Away  from  A u  ”0,  the  first  term  on  the  right- 
hand  side  of  Eq.  (4a)  is  the  only  contribution  to  the  total 
PDF 

An  alternative  expression  can  be  obtained  for 
Pto(Au(r))  as  follows.  Following  Kolmogorov  [13]  we 
assume  that  for  r  in  the  inertial  range  of  scales 

Au(r)-r(rc,),/3.  (5) 

where  r  is  a  “universal”  stochastic  variable,  and  re,  is  the 
total  energy  dissipation  in  the  linear  piece  of  size  r.  It  is 
clear  that  within  this  framework  [14]  any  reasonable 
model  for  e,  and  c  will  also  yield  a  reasonable  model  for 
the  PDF  of  A u(r).  A  convenient  model  [15]  for  the  ener¬ 
gy  dissipation  is  one  in  which  the  average  energy  flux 
summed  over  any  box  of  size  r/Lo~2~”  can  be  written 
as 


* 


rer-LoeoU  . 

/-i 


(6) 


AuOO-rfAu^^rim,1'3.  (7) 

l-l 

Recalling  again  from  Fig.  I  that  ^(Amo)  is  Gaussian 
(with  mean  0  and  variance  <Au^)'/2),  we  obtain  from  Eq. 
(7),  with  n“0,  that  the  universal  stochastic  variable  v 
has  a  normal  distribution  with  zero  mean  and  unity  vari¬ 
ance.  In  stating  this  result,  we  arc  stretching  the  validity 
of  Eq.  (5)  all  the  way  to  the  integral  scales  of  motion. 

Now,  Meneveau  and  Sreenivasan  [16]  have  shown  that 
the  following  simple  model  (the  so-called  p  model)  for 
the  probability  density  of  m  is  adequate  for  most  pur¬ 
poses: 

Pm(m)  —O.SlSlm  —  M)  +  5(m  —  (1  —  M)))  (8) 


with  M  ™0.3.  Using  the  normality  of  v  and  Eq.  (8),  we 
get  from  Eq.  (7)  that 


PsuiAufr))—  Z  Cf2 

*  -o 


(2*<rDl/2 


exp 


Am" 
2  <**;„ 
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(9a) 

(9b) 


where  the  multipliers  mt  are  identically  distributed  in¬ 
dependent  random  variables.  Here  eoLo  is  the  typical  en¬ 
ergy  dissipation  contained  in  an  eddy  at  the  macroscopic 
scale  (or,  equivalently,  the  total  energy  flux  across 
scales).  Noting  that  (Au$)~(Lo£o)Vi,  we  write 


Figure  2  shows  comparisons  between  the  experimental 
data  for  two  values  of  r  (r/A  “0.009  in  the  inertial  range 
and  r/A  “0.9)  and  the  theoretical  results  of  Eqs.  (4)  and 
(9).  The  appropriate  values  of  n  to  be  used  in  Eqs.  (4) 
and  (9)  are  determined  by  the  relation  n  “logj(Lo/r ). 
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FIG.  2.  A  comparison  between  the  experimental  data  for 
p*,(Au)  for  two  values  of  r,  and  Eq.  (9)  based  on  the  p  model 
(dashed  line)  and  Eq.  (4)  based  on  the  random  p  model  (dotted 
line).  The  formula  based  on  the  random  p  model  tends  to  un¬ 
derestimate  the  PDF  of  Au  close  to  0. 

Both  expressions  fit  the  data  reasonably  well,  but  the 
theoretical  curve  for  the  random  P  model  falls  below  the 
experimental  data  in  the  region  near  Au  “0.  This  is  so 
because  in  this  model  the  active  zones  (Au* 0)  tend  to 
become  more  sparse  as  n  increases. 

Equation  (9)  expresses  the  probability  of  finding  a  ve¬ 
locity  increment  Au  at  a  fixed  scale  r.  The  probability  of 
finding  a  “gradient”  s  “ Au/r  at  the  scale  r  can  be  writ¬ 
ten,  by  a  change  of  variables  from  Au  to  s,  as 

P(s\r)wmrpAu(rs)  ,  (10) 

where  pAu  is  given  by  Eq.  (9)  or  Eq.  (4).  The  quotation 
marks  for  the  gradient  above  reflect  the  fact  that  s  is  the 
gradient  only  when  r  is  comparable  to  the  dissipation 
scale  rD  defined  by 

roA«(ro)/v”I  .  (11) 

Note  that  the  average  value  (rD)  of  rD  is  the  Kolmogorov 
scale  rj.  For  the  p  model,  Eq.  (II)  becomes 

v((Aui)mLo/v)Mk,i(\ -M)(N~k)n 2_/v-l  .  (12) 

Here,  we  have  used  the  relation  rp/Lo~2~N,  and  k. 


FIG.  3.  A  comparison  between  the  experiment  for  the  tails  of 
the  p,(s)  (solid  line)  and  the  theoretical  formulas  based  on  the 
p  model  (dashed  line)  and  the  random  p  model  (dotted  line). 


FIG.  4.  The  stretching  exponent  m  determined  empirically 
for  Eqs.  (4)  and  (9),  compared  with  the  experimentally  deter¬ 
mined  data  of  Fig.  I  (dashed  line,  p  model;  dotted  line,  p  mod¬ 
el).  The  parameter  n  is  determined  arbitrarily  by  matching  the 
experiment  with  the  theoretical  formulas  for  some  r/A.  Inset: 
An  example  in  which  Eq.  (4)  for  n“IO  (corresponding  to 
r/A “0.009)  is  well  fitted  by  a  stretched  exponential  with 
stretching  exponent  m  “0.9. 


which  lies  between  0  and  N ,  is  randomly  picked  for  the  p 
model  from  a  binomial  distribution.  Equation  (12)  and 
the  knowledge  of  the  PDFs  of  r  and  k  will  allow  us  to 
compute  the  probability  that  the  dissipation  scale  is  rD. 
Let  P,D(r)  be  the  probability  that  the  dissipation  scale  is 
r.  The  probability  density  for  the  gradient  can  now  be 
written  from  Eq.  (10)  as 

P.t(s)mmJ drP(s\r)Pro(r) ,  (13) 

where  we  have  weighted  the  conditional  expectation 
P(s\r)  by  the  probability  that  the  dissipation  scale  is  r. 
This  step,  while  rigorous,  would  necessitate  cumbersome 
computations.  We  avoid  them  here  by  making  the  sim¬ 
plifying  assumption  that 

p,(s)~P(s\r—rD) .  (14) 


FIG.  5.  Demonstration  that  the  PDFs  of  increments  w>, 
which  b  the  velocity  u  from  which  the  low-frequency  com¬ 
ponents  are  removed,  tend  to  an  exponential  form  for  large  r. 
This  is  in  contrast  to  the  unfiltercd  data  for  which  the  asymp¬ 
totic  form  b  nearly  Gaussian.  The  high-pass-filter  setting  is  30 
Hz  for  +  and  9  Hz  for  0 .  The  error  bars  are  comparable  to 
the  typical  ones  shown  in  Fig.  I . 
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Noting  that  rp  ““  ( v/j ) ,/2,  Eqs.  (10)  and  (14)  allow  us  to  express  p,(s)  for  the  p  model  as 
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(15) 


(to  within  a  normalization  constant).  Note  that,  in  Eq. 

(15),  /V(s)  “  y  log2(Ltf|s|/v')  and  that  the  approximate 
sign  in  Eqs.  (14)  and  (15)  is  a  reminder  about  their 
heuristic  nature.  Using  these  same  assumptions,  Benzi  el 
al.  l9l  derived  a  similar  expression  for  the  PDF  of  the  ve¬ 
locity  gradient  on  the  basis  of  the  random  p  model.  Fig¬ 
ure  3  shows  a  comparison  between  experiment  and  Eq. 
(15)  and  an  alternative  formula  due  to  Benzi  el  al.  {9], 
The  p  model  yields  a  slightly  better  (it  to  the  data. 

We  now  point  out  that  the  stretched  exponentials  ap¬ 
proximate  Eqs.  (4)  and  (9)  even  though  the  latter  two 
are  substantially  more  complex.  This  is  illustrated  in  the 
inset  to  Fig.  4  where  Eq.  (9)  for  n“  10  is  plotted  against 
(Am)09.  The  tails  of  the  PDF  are  closely  approximated 
by  a  stretched  exponential,  consistent  with  Figs.  I  and  2. 
The  theoretical  expressions  have  been  examined  for  vari¬ 
ous  n,  and  the  stretching  exponents  m  are  determined 
from  the  best  empirical  (its  [17].  For  the  stretched  ex¬ 
ponentials  pAu(A«),“PA«(0)exp(  —  a|Au|'”),  reasonably 
good  straight  lines  can  be  fitted  to  the  plot  of 
log{log[p4„(0)//74l((Au)]}  vs  log(|Au|)  in  the  range  of 
|Au|  between  the  minimum  and  maximum  of  [see 
Eqs.  (4b)  and  (9b)].  The  exponents  so  determined  are 
plotted  in  Fig.  4  by  matching  with  experimental  data  at 
one  chosen  location.  The  fits  are  good  on  the  whole. 

We  now  return  to  the  asymmetry  of  the  PDF  of  A u(r). 
In  an  earlier  work  [7],  the  velocity  signal  u  was  decom¬ 
posed  into  its  Fourier  modes,  and  the  modes  below  the 
low-frequency  end  of  the  inertial  range  were  eliminated. 
The  remaining  Fourier  modes  were  recombined  to  give 
the  filtered  signal,  say  u  > .  The  increments  A u  >  (r) 
“u  >  (x  +r)  —  u  >  (x)  were  then  found  to  be  symmetri¬ 
cally  distributed.  By  fitting  stretched  exponentials  to  the 
tails  of  the  p^(Au  >  ),  we  found  that  the  exponents  were 
similar  to  those  of  the  unfiltered  data  in  the  dissipation 
range,  but  tended  towards  an  exponential  for  large  r  (Fig. 
5).  An  inquiry  on  the  exponential  behavior  of  the  tails  of 
the  PDFs  of  the  temperature  field  has  recently  been  un¬ 
dertaken  [18],  but  it  is  not  clear  whether  similar  explana¬ 
tions  hold  for  velocity  increments  in  the  filtered  signal. 
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